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Abstract 

We have measured the annual parallax of the H2O maser source associated with an infrared dark cloud 
MSXDC G034.43+00.24 from the observations with VERA (VLBI Exploration of Radio Astrometry). The 
parallax is 0.643±0.049 mas, corresponding to the distance of 1.56lg J 2 kpc. This value is less than the half 
of the previous kinematic distance of 3.7 kpc. We revise the core mass estimates of MSXDC G034. 43+00. 24, 
based on virial masses, LTE masses and dust masses and show that the core masses decrease from the 
previous estimations of ~ lOOOM© to hundreds of M@. The spectral type derived from the luminosity also 
changes from 09.5 to Bl in the case of MM1. This spectral type is still consistent with that of the massive 
star. The radial velocity derived from the flat rotation model is smaller than the observed velocity, which 
corresponds to the peculiar motion of ~ 40 km s -1 in the line-of-sight direction. 

Key words: astrometry — ISM: individual (MSXDC G034. 43+00. 24) — stars: formation — tech- 
niques: interferometric 
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1. Introduction 

Infrared dark clouds (hereafter, IRDCs) arc discov- 
ered by the observations with ISO (Infrared Space 
Observatory) and MSX (Midcoursc Space Experiment) 
(Perault et al. 1996; Carey et al. 1998; Hennebelle et al. 
2001), and observed as dark silhouettes against the back- 
ground radiation of our galaxy in the mid infrared wave- 
length (Perault et al. 1996; Egan et al. 1998). They are 
massive (~ 1O 2 -1O 4 M ), dense (> 10 5 cm -3 ), high col- 
umn density (~ 10 23 -10 25 cm -2 ), and low temperature 
(< 25 K) (Egan et al. 1998; Carey et al. 1998; Carey et al. 
2000; Tayssier, Hennebelle, Perault 2002; Simon et al. 
2006a; Simon et al. 2006b). 

Recently, IRDCs have been known as the sites where 
massive star formation is active (Rathborne, Jackson, 
Simon 2006; Rathborne, Simon, Jackson 2007; Pillai et al. 
2006; Jackson et al. 2008; Chambers et al. 2009). Various 
parameters, such as masses and luminosities, have been 
obtained from many observational studies from the mid- 
infrared to millimeter wavelengths (e.g., Rathborne et al. 
2005; Shepherd et al. 2007; Sanhueza et al. 2010). These 
parameters depend on distance, and in most cases, kine- 
matic distances have been used to derive the parameters 
assuming the rotational model of our Galaxy. However, 
recent annual parallax measurements with VLBI show 
that kinematic distances tend to be overestimated and 
sometimes more than the double of those from parallaxes 
(e.g. Sato et al. 2010a, Motogi et al. 2011). Therefore, 



it is important to measure distances toward the sources 
accurately to know the physical condition of the sources. 

Annual parallaxes of galactic sources located at a few 
kpc from the Sun can be measured from the phase- 
referencing VLBI observations. Phase-referencing VLBI 
is technique to observe two adjacent (less than several de- 
grees in the case of 22 GHz-band observations) sources 
on the sky simultaneously or in the short time interval 
(< 60 s at 22 GHz band) (Asaki et al. 2007; Thompson, 
Moran, Swenson 2001). By subtracting the difference of 
optical path lengths between the the target and reference 
source, it becomes possible to eliminate the short-term 
fluctuations of the earth's atmosphere. In addition, the 
reference position is given by selecting an extragalactic 
source as the reference source. These benefits enable the 
astrometry referred to the reference source and, with the 
long-time integration, the detections of faint sources which 
is not able to detect with the normal VLBI observations. 

Phase-referencing VLBI observations with VERA 
(VLBI Exploration of Radio Astrometry) and the VLBA 
(Very Large Baseline Array) derive annual parallaxes of 
the kilo-parsec scale for the galactic maser sources (e.g., 
Baba et al. 2009 and references therein; Sato et al. 2010a; 
Oh et al. 2010; Rygl et al. 2010; Kamohara et al. 2010) af- 
ter the first successful measurements of annual parallaxes 
with H2O mascrs by Kurayama, Sasao, Kobayashi (2005). 
It is possible to measure the annual parallaxes when both 
the maser sources and the adjacent extragalactic reference 
sources are bright and compact sufficient to observe with 
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VLBI. 

Some IRDCs emit masers such as the H2O masers 
and CH3OH masers (Pillai ct al. 2006; Wang ct al. 
2006; Chambers ct al. 2009), so wc can measure the 
annual parallaxes from the phase-referencing VLBI ob- 
servations of these masers. We used H 2 masers this 
time because most parallaxes measured with VERA are 
based on the H2O maser observations. Wang et al. 
(2006) surveyed H2O masers for IRDCs and found that 
MSXDC G034.43+00.24 has bright H 2 masers. There 
is also an extragalactic source as a reference source for 
phase-referencing within 2° from MSXDC G034.43+00.24 
(Fomalont et al. 2003). Therefore, this source is one of 
the most suitable IRDCs to measure the parallaxes with 
VERA. 

MSXDC G034. 43+00. 24 has a filamentary structure ex- 
tended over ~ 9' in the declination direction, and four 
millimeter cores (MM1-MM4) on this filamentary struc- 
ture (Faundcz et al. 2004; Garay et al. 2004; Rathborne 
et al. 2005). An IRAS point source (IRAS 18507+0121) 
and an ultra compact H II region are associated with MM2 
(Miralles, Rodrigues, Scalise 1994; Molinari et al. 1998). 
MM1, MM3 and MM4 have 4.5 (im excesses, which sug- 
gest the existence of the ionized or shocked gas. Point 
sources seen at 24 fim and associated with all four mil- 
limeter cores suggest the existence of the warm dust 
(Rathborne et al. 2005; Chambers et al. 2009). Garay 
et al. (2004), Rathborne et al. (2005), Rathborne, Jackson, 
Simon (2006), Shepherd et al. (2007) and Sanhueza ct al. 
(2010) calculated the mass of each millimeter core. VLA 
observations have detected H2O masers in MM1, MM3 
and MM4 (Wang et al. 2006). No VLBI observation has 
been conducted for MSXDC G034.43+00.24. The avail- 
able distances to MSXDC G034.43+00.24 is the kinematic 
distance only, whose value is 3.7 kpc, from 13 CO J = 1 — > 
observations (Simon et al. 2006b) and CS(2-1) observa- 
tions (Bronfman, Nyman, May 1996; Faundez et al. 2004). 
Hence, it is very important to measure the annual parallax 
of this source. Here, we report the results from VERA. 

2. Observations 

Observations were carried out with VERA four sta- 
tions. The target maser is the 616 — > 5 2 3 transition of 
H 2 molecules, whose rest frequency is 22.23508 GHz. 
We monitored at ten epochs from 2006 November to 2008 
July with the typical observing duration of 9 hours. The 
observation dates are listed in table 1. Before this astro- 
metric monitoring, we conducted maser survey observa- 
tions of four millimeter cores on 2006 October 21. 

We conducted phase-referencing VLBI observations of 
the target source MSXDC G034.43+00.24 and the refer- 
ence source GPSR5 35.946+0.379 (= VCS2 J1855+0251, 
hereafter we call this source as the reference source), 
which is separated by l.°6 from the target source. Using 
the dual beam system of VERA telescopes (Kawaguchi, 
Sasao, Manabe 2000; Honma ct al. 2008), wc observed 
these two sources simultaneously. For the fringe finding 
and bandpass calibration, we also observed the calibra- 



tor source QSO J1800+3848 (a J2 ooo = 18 h 00 m 24?765362, 
<$J20oo = +38°48'30"69754) with the both beams every 80 
minutes. 

Observed signals are filtered with the digital filter bank 
(Iguchi et al. 2005) and outputted in the form of six- 
teen IFs. Each IF contains 16-MHz bandwidth signal. 
One IF contains the signal of the target (maser) source. 
The other fifteen IFs contain the signals of the reference 
(continuum) source. The recording rate is 1024 Mbps 
(32 Msps x 2 bit sample -1 x 16 IF). The correlation is 
performed with Mitaka FX correlator. The correlator in- 
tegrates the output data for 1 s in the time domain. In 
the frequency domain, the output data are 8 MHz band- 
width with the 512 frequency points for one maser IF and 
16 MHz bandwidth with the 64 frequency points for four- 
teen continuum IFs. One continuum IF is abandoned in 
the data reduction process. The velocity resolution at 
the correlator output is 0.21 km s _1 for the maser IF and 
3.4 km s" 1 for the continuum IFs. 

The measurements of the annual parallax of MSXDC 
G034. 43+00. 24 require a special care because this source 
has a low-declination of ~ +1°. This causes a decrease 
of the positional accuracy in the declination direction. In 
VLBI observations, we vary the measurement points on 
the uv plane, which corresponds to the aperture of the nor- 
mal single-dish telescopes, by changing the positional re- 
lationship between the stations and the sources using the 
earth's spin. The loci of measured points on the uv plane 
draw parts of the ellipses for high-declination sources, but 
they draw almost straight lines in the right-ascension di- 
rection for low-declination sources. Thus the positional 
accuracy in the declination direction decreases especially 
in the case of the small number of stations, that is, small 
numbers of baselines and observing points in the decli- 
nation direction. In our parallax measurement, we used 
only the right-ascension components in order to avoid this 
effect. The reason of large uncertainty in the declination 
direction is also described in §4.3.1. 

3. Data Reduction 

Data reduction was done with NRAO (National Radio 
Astronomy Observatory) AIPS (Astronomical Image 
Processing System) software and DIFMAP. The functions 
neither AIPS nor DIFMAP provide were added with ex- 
ternal programs. The steps of our data reduction are as 
follows : 

1. Amplitude calibration for all sources 

2. Bandpass calibration for all sources using calibrator 
source 

3. Modification of tracking model in Mitaka FX corre- 
lator 

4. Time and frequency integration 

5. Calibration of clock parameters (Global fringe 
search for calibrator source) 

6. Global fringe search for reference source 

7. Imaging of reference source with self-calibration 

8. Phase referencing (Subtracting the phases of refer- 
ence source from those of target source) 
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Table 1. Observation summary. 
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epoch 


date 


day offset* 


synthesized beanr 


1st 


2006 Nov 16 


-297 


1.29 x 0.86 mas at 


-53° 


2nd 


2007 Jan 04 


-248 


1.40 x 0.79 mas at 


-50° 


3rd 


2007 Feb 10 


-211 


1.39 x 0.82 mas at 


-52° 


4th 


2007 Mar 25 


-168 


1.36 x 0.81 mas at 


-52° 


5th 


2007 May 04 


-128 


1.28 x 0.74 mas at 


-53° 


6th 


2007 Aug 06 


-34 


1.32 x 0.81 mas at 


-44° 


7th 


2007 Oct 02 


23 


1.33 x 0.80 mas at 


-45° 


8th 


2008 Jan 04 


117 


1.43 x 0.79 mas at 


-49° 


9th 


2008 Apr 12 


216 


1.35 x 0.82 mas at 


-48° 


10th 


2008 Jul 02 


297 


1.27 x 0.81 mas at 


-39° 



* Offset from the center date (2007 September 9) of the term. 

t The lengths of major and minor axes and the position angle of the major axis. 



9. Dual-beam calibration with 'horn-on-dish' method 
(Honma et al. 2008) 

10. Calibration of Doppler effect 

11. Imaging of target source with CLEAN 

12. Measuring the position and flux by fitting elliptical 
Gaussians 

The details of the data reduction procedure for an observ- 
ing epoch are shown in Appendix 1. 

Modification of tracking model in Mitaka FX cor- 
relator' (step 3) is needed because the tracking model 
in Mitaka FX correlator does not have sufficient accu- 
racy for astrometric measurements. We changed it to 
the newly calculated model with CALC3/MSOLV (Jike 
et al. 2005; Manabe, Yokoyama, Sakai 1991). In this 
recalculation, we used the equatorial coordinates for 
the target source MSXDC G034.43+00.24 as aj 20 oo = 
18 h 53 n T9?00000, c5j2ooo = +1°24'08."0000 and for the ref- 
erence source GPSR5 35.946+0.379 = VCS2 J1855+0251 
as aj2ooo = 18 h 55 m 35. s 43649, <5j 20 oo = +2°51'19."5623 
(Fomalont et al. 2003). The equatorial coordinate for 
the target source corresponds to that of the millimeter 
core MM1 (Rathborne et al. 2005). The new recalculated 
model contains the estimation of the wet component of 
the earth's atmosphere from GPS data at VERA stations 
(Honma, Tamura, Reid 2008). 

Before 'imaging of target source with CLEAN' (step 
11), we need to know the position where we make images. 
Only the H2O masers associated with the millimeter core 
MM1 was imaged (see §4.1). Producing the fringe-rate 
map with our own program from the results of the global 
fringe search for the flux-peak frequency chennel, we find 
a maser cluster as shown in figure 1. We made the im- 
ages for the eight features in the maser cluster at which 
we found the masers for the data on 2008 January 4. We 
adopted maser features only when their radiation is con- 
firmed over three or more frequency channels. 

In 'measuring the position and flux by fitting elliptical 
Gaussians' (step 12), we fit the elliptical Gaussians for the 
five regions where we detected maser emission during more 
than two continuous epochs (see, table 2). We adopted the 
peak positions as the positions of maser features. In this 
paper, maser features denotes maser emissions detected 




RA offset [mas] 

Fig. 1. Fringe rate map of the H2O masers at the mil- 
limeter core MM1 made from the results of the global 
fringe search with AIPS. The map origin is located at 
the position of the target source modified in the recalcula- 
tion of the tracking model whose equatorial coordinate is 
QJ2000 = 18 h 53 m 19f 00000, <5 J20 oo = +1°24'08"0000. Each 
line shows the possible position of the masers calculated from 
one rate value. Maser clusters are placed at the crossing 
point of lines. This map is obtained from the data on of 
2008 January 4 at the radial velocity of ulsr = 56.5 km s _1 . 

at multiple frequency points at the same position on the 
sky. Positions are measured at the frequency channels at 
which the flux is largest among the detected frequency 
channels. Detected maser features are compact sufficient 
for the position measurements, shown in Figure 2. When 
multiple maser features are detected in an image at an 
epoch, we used the feature which has closest declination 
value to the other epochs. Figure 3 shows an example of 
the selection of maser features in this process. From this 
figure, we can find that declinations are better for this 
selection than right ascensions although declinations have 
large position errors shown in §4.3.1. 
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Fig. 4. Cross-power spectra of the H2O masers at the millimeter-core positions: (a) MM1, (b) MM2, (c) MM3, and (d) MM4. 
They are obtained from the data on 2006 October 21 at Mizusawa-Iriki baseline. The integration time of data is 5 minutes. We 
integrated the data after the phase-referencing, which gets rid of the short-time fluctuation of the earth's atmosphere and enables 
the long-term integration. 



4. Results 

4-1. Survey Results of Masers 

Figure 4 shows the cross-power spectra obtained from 
the survey observations on 2006 October 21 at the four 
millimeter-core positions. We find that the H2O masers 
are detected in MM1 and MM4, but are not detected in 
MM2 and MM3. Comparing with the VLA observations 
in 2006 (Wang et al. 2006), the detection is not consistent 
only for MM3. This may be caused by the time variation 
of masers or the extended emission of the source structure, 
that is, masers in MM3 could be resolved out with the long 
(~ 2000 km) baselines of VERA. 

From this result, we observed the millimeter-core MM1 
for the annual parallax measurements from the observa- 
tion on 2006 November 16. 

4-2. Derivation of Parallax with Least Square Fitting 

Figure 5 and 6 show the movements of maser features 
measured with respect to the reference source. The 
movements look like the combination of the parallax mo- 
tion and the linear motion in the right-ascension direction, 
but have large scatter over 20 mas in the declination direc- 
tion. This is because the declination of the target source 
(~ +1°) is small. As shown in figure 7, there are large 
sidelobes of the synthesized beam in the declination di- 
rection. The stable and uncalibrated component of the 



earth's atmosphere affect the positions in the declination 
direction. The details of this effect is discussed in §4.3.1. 
We used only the data of the right-ascension direction this 
time. 

We carried out the least square fitting only with the 
right-ascension data to derive the annual parallax. The 
model is the combination of the movement by the annual 
parallax and the linear motion. The movement except 
by the annual parallax, such as the galactic rotation, the 
peculiar motion of the Sun, arc contained in the linear 
motion. This is expressed as the following equation: 

or*' (t) cos 5 = w(— sinacos ® + cosecosasin ® ) 



+ (n® cos5)t + a§' cos 8 



(1) 



where cr^(i) is the observed right ascension of the i'th 
maser feature at the time t, a and 6 are the right ascension 
and declination of the target source respectively, w is the 
annual parallax, ® is the ecliptic longitude of the Sun and 
e is the obliquity of ecliptic. fj,a coaS is the right-ascension 
component of the linear motion of the i'th maser feature, 
t is the day offset from the center of the observation term 
used for the parallax measurements, shown in table 1 and 
(Xq is the right ascension of i'th maser feature at the time 
t = (Green 1985). 

Table 2 summarizes the detection of masers at each fea- 
ture and at each epoch. Maser features detected during 
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Table 2. Detection at each feature and epoch. 
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^LSR 


RA* 


dec* 






Integrated flux at each 


epoch [Jy] T 






Number 


[km s 1 


[arcsec] 


[arcsec] 


1st 


2nd 


3rd 


4th 


5th 


6th 


7th 


8th 


9th 


10th 


i 


88.7 


-3.632 


+0.088 




(1.4) 


(1.1) 








(2.6) 


(3.6) 






2 


61.7 


-3.489 


+0.819 


(3.4) 


5.9 


8.8 


3.8 


(5.7) 


3.1 


2.5 


9.4 


(2.8) 




3 


59.4 


(-3.367) 


(-0.036) 




(2.1) 










(5.2) 




(2.2) 




4 


57.9 


-3.425 


-0.061 


(2.4) 


(4.4) 


(3.9) 








(0.8) 


(2.0) 






5 


57.3 


-3.303 


-0.035 


(1.7) 


(3.0) 










1.0 


(5.3) 


22.4 


11.5 


6 


56.5 


-3.672 


+0.107 


(3.2) 


(6.6) 


(5.8) 


(4.5) 






(2.6) 


(27.3) 


(54.5) 




7 


53.9 


-3.425 


+0.657 








2.2 


(8.7) 


3.8 


8.7 


19.2 


(5.8) 




8 


43.4 


-3.392 


-0.048 
















(5.1) 


(0.5) 





* Offset values from the right ascension and declination used in the modification of the tracking model : aj2000 = 18 53 m 19?0, <5j2000 = 
+1°24'08". Values without parentheses are the offsets at 8th epoch on 2008 January 4. Values with parentheses are the averages of the 
detected epochs. 

t Numbers without parentheses denote that the data of the features and epochs are used for the annual parallax measurements. Numbers 
with parentheses denote that the data are not used. 
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Fig. 2. An example of the phase- referenced images of de- 
tected maser features. The phase-referenced image of maser 
feature No. 2 at 3rd epoch on 2007 Feburuary 10. Contours 
are 3cr, 5cr, la and 9rr brightness level, where a is the rms 
noise level of 0.99 Jy beam" 1 . Maser features are compact 
sufficient to measure their positions. 
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Fig. 3. Maser-feature selection for the feature No. 2 when 
multiple maser features are detected in an image at an epoch, 
(a) Plot of right ascensions versus time, (b) Plot of declina- 
tions versus time. In both panels, open circles (O) show the 
maser features which are used for the parallax fitting, and 
crosses (x) show the maser features which are not used for 
the parallax fitting. From (a), it is difficult to know which 
feature should be used, but we can easily select from (b) al- 
though declination data have large position errors shown in 
§4.3.1. 
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Fig. 5. Observed data points in the right ascension direction 
and the least square fitting for the annual parallax measure- 
ments. Table 1 shows the correspondence between the values 
of day offsets from 2007 September 9 (horizontal axis) and 
the observation dates. The offset in the vertical axis for each 
feature is included just for the display. The observed points 
are the center of the error bars. Crosses (x) show the data 
points which are not used for the parallax fitting. The ampli- 
tudes of uncertainty are O.lOf mas, which is calculated from 
the assumption of (reduced x 2 ) = X 2 l u — 1- Feature numbers 
are same as those in table 2. Solid curves are the results of the 
least-square fitting as shown in table 3. Dashed curves show 
the movements of features which are derived from the least 
square fitting for calculating linear motions and initial posi- 
tions only, assuming the parallax of vj = 0.643 mas derived 
from the fitting. 



* 
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Fig. 7. The synthesized beam of VERA four stations in the 
observation on 2008 January 4. Strong sidclobes are found in 
the declination direction. 



Table 3. Results of the fitting of the annual parallax and the 
linear motions 




W 



0.643 ±0.049 mas 



-300 -200 -100 100 200 
Day offset from 2007 Sep 09 [days] 



300 



Fig. 6. Plot of the observed data points in the declination 
direction against time. Table 1 shows the correspondence be- 
tween the day offsets from 2007 September 9 (horizontal axis) 
and the observation dates. The offset in the vertical axis is 
included for each feature for the display. Crosses (+) show the 
positions of feature number 2, asterisks (*) of feature num- 
ber 4, circles (O) of feature number 5, Crosses (x) of feature 
number 6 and squares (□) of feature number 7. Feature num- 
bers are same those in table 2. Solid, dotted and dashed 
curves show the movements of feature number 2, 7 and 5 re- 
spectively, which are derived from the least square fitting for 
calculating linear motions and initial positions only, assuming 
the parallax derived from right-ascention data vj = 0.643 mas. 



Feature^ 
Number 



IT 

jJL a COS0 + 

[mas yr -1 ] 



1*7 



'o 



COS(0 



mas 



-3488.77 ±0.05 
-3303.80 ±0.12 
-3425.55 ±0.05 



2 -2.19±0.12 
5 0.98 ±0.21 

7 0.46 ±0.19 

* Value ot the annual parallax. 

* Corresponding to the feature numbers in table 2. 

* Right-ascension components of the linear motions. 

§ Offsets in the right-ascension direction at t = on 2007 
September 9 from the position used in the modification of tracking 
model ajaooo = 18 h 53 m 19?0, <5 J2 ooo = +1°24'08". 



more than two continuous epochs are used for the par- 
allax measurements. The data points which are not on 
the declination trend are omitted from the fitting. This 
is because we can not solve parameters when we can not 
use data of more than two epochs. When we include one 
new maser feature, we need two extra fitting parameters: 
the initial position and the linear motion. We fitted for 
three maser features, so seven parameters are fitted : one 
common w, five /a« , five q;q^cos<5. The total number of 
data points is twenty-nine, as shown in table 2. 

Figure 5 and table 3 show the results of fitting. The de- 
rived annual parallax is 0.643 ±0.049 mas, corresponding 
to the distance of 1.56^°'^ kpc. This distance is smaller 
than the half of the previous kinematic distance of 3.7 kpc. 
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where B(x, y) is the brightness distribution of the ref- 
erence source, x and y are the coordinates of the right- 
ascension and declination direction respectively, u> = 2ttu 
and v is the observing frequency (Fujishita 1983; Chariot 
1990). For the simple calculations, we set B(x,y) as 
the sum of the point sources obtained from CLEAN, 
B(x,y) ^J2i^oiS(xi,yi), and approximate d/du:^ 1/cj = 
1/(2tti/) and 2ttut(u, v, cj) ~ 2tt\/u 2 + v 2 Ax, where Ax 
is the positional error. Then, so we approximate 
27tz/t(m,w,cj) ~ 2-k\Ju 2 + v 2 Ax where Ax is the positional 
error. Then 



Ax 



o 



10 5 -5 

MilliARC SEC 
Center at RA 18 55 35.4364901 DEC 02 51 19.562320 
Cont peak flux = 1.4208E-01 JY/BEAM 
Levs = 1.421E-02 * (0.449, 0.897, 1.794, 3.588, 
7.176) 



Fig. 8. Imaging result of the reference source GPSR5 
35.946+0.379 = VCS2 J1855+0251 with VERA four stations 
at 22 GHz band. The observation date is 2008 January 4. 
Contours are 3cr, 6a, 12a, 24cr and 48<r brightness level, where 
a is the rms noise level of 2.1 mjy beam . 

4--S. Uncertainty of Each Epoch 

Uncertainty at one epoch is calculated by assuming that 
the uncertainty values of all epochs are same and that 
(reduced y 2 ) = \ 2 l v = L where \ 2 is the squared sum of 
the residuals of fitting over the square of the uncertainty 
of one epoch, v is the degree of freedom in the fitting, the 
total number of data points minus the total number of 
parameters. In our fitting, v = 29 — 11 = 18. The uncer- 
tainty at one epoch is 0.101 mas. The uncertainty of fitted 
parameters is calculated from this value with the normal 
least square fitting method. Most of the uncertainty is 
the stable component of the earth's atmosphere of about 
0.2 mas, as discussed in §4.3.1. 

The uncertainty from the structure of the reference 
source is small. Figure 8 shows the imaging result of the 
reference source. This source structure has little change 
over the all observation epochs used for the parallax mea- 
surements. We used same stations in all observing epochs, 
so the measurement points on the uv plane are basically 
same. 

In order to evaluate the uncertainty of the parallax mea- 
surements from the structure of the reference source, we 
calculated the worst value of the uncertainty as follows. 
The delay time, optical path length divided by the speed 
of light, caused by the source structure t(u,v,oj) is given 

by 



-(u,v,ui) = 



arg 



B{x lV )e- 2 ™ {ux+vy) dxdy 



■ (2) 



2tt\/^2" 



: arg 



^Boi 



-2iri(uxi-\-vyi) 



(3) 



We calculated this value for the data on 2008 January 4 
and the maximum values of u and v of w max = 1 .4 x 10 8 and 
fmax = 1-2 x 10 8 in the unit of the wavelength. The result 
is Ax = 0.045 mas. The uncertainty from the structure 
of the reference source has a small fraction in the each 
epoch's uncertainty of 0.101 mas. 

4-3.1. Uncertainty from Stable Component of the Earth's 
Atmosphere 

In the phase-referencing VLBI observations, the effect 
of the earth's atmosphere is removed by taking the differ- 
ences of the optical path lengths between the target and 
reference sources. With this process, we can remove the 
short-term fluctuations of the earth's atmosphere, which 
is the largest atmospheric effect in VLBI observations. 
However, the elevation angle of the target and reference 
sources are not completely same, so it is impossible to re- 
move the effect of airmass, that is, the stable component 
of the earth's atmosphere (Honma, Tamura, Reid 2008). 
Assuming the plane-parallel approximation of the earth's 
atmosphere, the effect of airmass is given by in the unit 
of the optical path length 

lij Lj sec z^j , (4) 

when source i is observed with the station j, where Lj 
is the airmass effect in the zenith direction in the unit 
of the optical path length, and Zij is the zenith angle of 
the observed source. In phase-referencing observations, 
we subtract the optical path lengths between the target 
and reference sources. After the phase-referencing, the 
following airmass effect remains: 



hj — hj = Lj (sec zij — secz2j) 



(5) 



We evaluated quantitatively how much this effect con- 
tributed to the uncertainty of our astrometric observations 
in the following way: (1) Giving the optical-path-lcngth 
offset in the zenith direction by the earth's atmosphere, L. 
We added this offset to Ishigakijima station only for the 
simplicity because this station is most humid and hence, 
L could have large uncertainty. We changed L at 1 cm 
step in the range of —30 cm < L < 30 cm. (2) Adding the 
phase correction of 

2ir 

— L[sec z t (t) — sec z r (t)] (6) 
A 
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obtained from the plane-parallel approximation to all vis- 
ibilities of the target source at the baselines including 
Ishigakijima station, where A is the observation wave- 
length, z t (t) and z r (t) are the zenith angle of the target 
and reference source at the time t, respectively. We used 
the relationship between the phase <p and the optical path 
length I of <j)= (2tt/X)1. (3) Imaging the target source from 
the phase-corrected visibilities without CLEAN, and de- 
riving the flux and position of the maser feature by picking 
up the pixel of the flux peak. 

Figure 9 shows the result of above evaluation for the 
maser feature in the feature number 5 and for the ob- 
servation on 2008 January 4. As discussed by Honma, 
Tamura, Reid (2008), the true atmospheric path length 
should maximize the peak brightness. There are not so 
large offset in the right-ascension and declination direc- 
tion at the atmospheric path length of about —8 cm, but 
the offset becomes large, more than 10 mas in the decli- 
nation direction at the atmospheric path length of around 
—22 cm. Such a large L would not be expected if GPS 
calibration is included (Honma, Tamura, Reid 2008) in 
which case it is normally a few centimeters. Therefore, 
the stable component of the earth's atmosphere might be 
one of the largest causes of the large position offsets in the 
declination direction seen in figure 5, but we can not deny 
other causes, for example the structure effect of the target 
maser source, large offsets (~ 4") from the right ascension 
of the tracking model. This atmospheric effect is the most 
part of the uncertainty in right ascension direction at each 
epoch of 0.101 mas, because the offset value in the right- 
ascension direction in the range of —25 cm < L < cm is 
about ±0.1 mas. 



4-4- Other Approaches for 
Measurements 



Poor Declination 



From figure 9, we can find that right ascensions are 
also affected by the stable component of the earth's at- 
mosphere. In order to evaluate this, we carried out the 
fitting of the annual parallax with the following six meth- 
ods: 

1. Using five features (feature numbers 2, 4, 5, 6 and 
7). Using all epochs, also including the data points 
which are far from the declination trends (1st, 5th 
and 9th epochs for feature number 2; 8th epoch for 
feature number 5; 5th and 9th epochs for feature 
number 7) except the 1st and 2nd epochs of feature 
number 5 which have a long time gap against the 
later epochs. Using the data of right ascension only. 

2. Using three features which have trends in the dec- 
lination measurements (feature numbers 2, 5 and 
7). Using all epochs, also including the data points 
which are far from the declination trends (1st, 5th 
and 9th epochs for feature number 2; 8th epoch for 
feature number 5; 5th and 9th epochs for feature 
number 7) except the 1st and 2nd epochs of feature 
number 5 which have a long time gap against the 
later epochs. Using the data of right ascension only. 




-30 -20 -10 10 20 30 
zenith path length at Ishigakijima station [cm] 



Fig. 9. (a) Right-ascension offset, (b) declination offset, (c) 
brightness of the peak-flux pixel against the optical path 
length of the zenith direction by the earth's atmosphere. In 
(a) and (b), the right ascension and declination not added the 
atmospheric effect arc set to zero. We evaluated for the maser 
feature in the feature number 5 observed on 2008 January 4. 
Since the pixel size is 0.1 mas and since the resolution of 
right-ascension and declination is the pixel size, the quan- 
tized variation in (a) is caused by the pixel size. The nature 
of quantized variation in (b) may be the effect of the sidclobc 
pattern (sec figure 7). 
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3. Using three features which have trends in the decli- 
nation measurements (feature numbers 2, 5 and 7). 
Excluding the epochs at which the declination data 
are far from the trends (1st, 5th and 9th epochs for 
feature number 2; 8th epoch for feature number 5; 
5th and 9th epochs for feature number 7) and the 
1st and 2nd epochs of feature number 5 which have 
a long time gap against the later epochs. Using the 
data of right ascension only. This is the method we 
adopt at §4.2. 

4. Same as the method No. 3, but using the data of 
both right ascention and declination. 

5. Using three features which have trends in the dec- 
lination measurements (feature numbers 2, 5 and 
7). Using all epochs, also including the data points 
which are far from the declination trends (1st, 5th 
and 9th epochs for feature number 2; 8th epoch for 
feature number 5; 5th and 9th epochs for feature 
number 7) except the 1st and 2nd epochs of feature 
number 5 which have a long time gap against the 
later epochs. For the epochs at which the declina- 
tion data are far from the trends, we tried to recover 
the declination by setting the CLEAN boxes at the 
likely sidelobcs from the trends of declination data. 
In two epochs, 5th and 9th epochs for feature num- 
ber 7, we failed this recovery, so we omitted them 
from the fitting. Using the data of right ascension 
only. 

6. Same as the method No. 5, but using the data of 
both right ascention and declination. 

Table 4 summarizes the result of these fitting. Basically, 
all fitting results of annual parallax coincide in the error 
range. Because the method No. 3 has smaller error than 
the methods No. 1, 2 and 5, right-ascension data also have 
large errors when their declination data have large errors. 
Comparing the erros of methods No. 3 and 5, we find that 
we can not recover these large-declination-error data by 
setting the CLEAN boxes at the likely sidelobes from the 
trends on declination data. Comparing method No. 3 and 

4, or method No. 5 and 6 show that the declination data 
have larger errors than the right ascension data even when 
we try to recover declinations with the settings of CLEAN 
boxes. In our observations, it might be better not to use 
the declination data for the fitting of annual parallax. We 
adopted the fitting method No. 3 because it is the best fit 
and does not include bad data. 

5. Discussion 

5.1. Modification of Physical Parameters of Star 
Forming Regions 

The distance to the target source MSXDC 
G034. 43+00. 24 becomes less than half of the kine- 
matic distance of 3.7 kpc. This changes the parameters 
of the star forming regions in this infrared dark cloud. 
We discuss some examples of them. 

Sanhueza et al. (2010) estimates the virial mass of each 
core from the following equation, assuming the spherical 



shape of each core and ignoring the magnetic field and 
external forces: 



M v 



Mr. 



= B — 



R \ ( Av 



pc J \ km s 1 



(7) 



where Av is the mean line width, R is the radius of the 
cloud, B is the constant depending on the density profile 
of the cloud. They calculated with B = 210, corresponding 
to the uniform density. Virial mass M v - lr is proportional to 
the distance, because M v - U - is proportional to the radius. 

They also calculated the LTE (local thermody- 
namic equilibrium) masses from the observations of the 
13 CO(3 -> 2) and C ls O(3 -> 2) lines. Assuming that 
C 18 0(3 — > 2) is emitted under the LTE condition and this 
emission is optically thin, LTE mass -Mlte is calculated 
from the following formula: 



Mlte 



0.565 



(T e; 



;)( 



[H 2 /C ia O] 
2.72m H y V 3 - 8x 106 
0.88)exp(15.81/T ox ) 



D 
kpc 

l-ex P (-15.81/r ox ) JJ^ d ^ W 

where [i m is the mean molecular weight per H 2 molecule, 
hih is the mass of the hydrogen atom, [H2/C 18 0] is the 
abundance ratio of C 18 relative to H2, D is the distance, 
T ex is the excitation temperature, which is employed to be 
T cx = 30 K, Tis is the optical depth of the C 18 0(3-s-2) line, 
v is the radial velocity in the unit of km s _1 , f2 is the solid 
angle in the unit of arcmin 2 . They derived Tig from the 
ratio of the observed brightness temperatures of C 18 0(3 — > 
2) to 13 CO(3 -> 2), assuming [ 13 CO/C 18 0] = 7.6. LTE 
mass is proportional to the square of the distance. 

Rathborne, Jackson, Simon (2006) estimated the dust 
mass from the 1.2 mm continuum observations with the 
following equation: 



M, 



F V D 2 



dust — D !rr \i (9) 

where Md us t is dust mass, F v is the observed flux density, 
D is the distance, K v is the dust opacity per unit mass, 
-B„(Tdust) is the Planck function at the dust temperature 
Tdust- They adopted Ki.2mm = 1-0 cm 2 g _1 , gas to dust 
mass ratio of 100, and the dust temperature of Td us t = 
15 K. Dust mass is also proportional to the square of the 
distance. 

Table 5 shows the masses calculated from the previous 
distance D = 3.7 kpc and our distance D = 1.56 kpc. The 
resultant masses become one order small, and the most 
massive core is no more than thousand Mq. 

Rathborne et al. (2005) derived the bolometric lumi- 
nosities from the standard graybody fitting of the spectral 
energy distribution (SED) between the mid infrared and 
millimeter wave. They assumed that the cores are isother- 
mal and that their radii are 15". The parameters derived 
from the fitting are the emissivity index /3, the optical 
depth at 250 (im T25o Aim , the dust temperature Td us t, and 
the bolometric luminosity Lboi- Furthermore, the lumi- 
nosity is constant for the high-mass protostar during the 
evolution from the protostar to the main-sequence, so the 
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Table 4. Result of various fittings for the annual parallax 



No. 


Method 


RA/dec 


w* 


(j ' 


N* 


n§ 


J 








[mas] 


[mas] 








1 


five features, all epochs 


RA only 


0.714 ±0.060 


0.212 


29 


11 


18 


2 


three features, all epochs 


RA only 


0.624 ±0.060 


0.182 


19 


7 


12 


3 


three features, epochs on dec trend 


RA only 


0.643 ±0.049 


0.101 


13 


7 


6 


4 


three features, epochs on dec trend 


RA and dec 


0.553 ±0.083 


0.181 


26 


13 


13 


5 


three features, all epochs, forced CLEAN boxes 


RA only 


0.726 ±0.066 


0.177 


17 


7 


10 


6 


three features, all epochs, forced CLEAN boxes 


RA and dec 


0.625 ±0.095 


0.269 


34 


13 


21 



* Annual parallax 

T Uncertainty in one epoch 

t Total number of data points 

§ Total number of calculated parameters 



Degree of freedom in the fitting 



Table 5. Physical parameters of the star forming region MSXDC G034. 43+00. 24 



Millimeter Core 



Values at the Kinematic 
Distance of D = 3.7 kpc 
MM1 MM2 MM3 MM4 



Values at Our Distance 
of D= 1.56 kpc 
MM1 MM2 MM3 MM4 



Reference 



Virial Mass 
LTE Mass 
Dust Mass 



1130 1510 
330 1460 
1187 1284 



1370 



301 



253 



476 
59 
211 



637 
260 
228 



578 



54 



(1) 
(1) 

45 (2) 



w 

(3) 



32000 ••• 9000 
09.5 ■■■ B0. 5 

't he unit ot masses and luminosities are Mq and L( 



^bol 

Spectral Type 



12000 
B0.5 



5700 
Bl 



1600 
B3 



2100 
B2 



0, respectively. 

(2) Rathborne, Jackson, Simon (2006), (3) Rathborne et al. (2005) 



Reference : (1) Sanhueza et al. (2010), 



bolomctric luminosities tell us the corresponding spectral 
types of protostars. These values are also shown in table 
5. They are still early B type star with the mass ~ 10M Q . 

5.2. The Position and Motion on the Galactic Plane 

MSXDC G034.43+00.24 is located at the closer part 
of Sagittarius-Carina arm. Figure 10 shows the position 
of MSXDC G034.43+00.24 at the distance of 1.4 kpc de- 
rived from our annual parallax. The background is cited 
from Sato et al. (2010a). The black lines and points are 
the model of the Galaxy by Georgelin, Georgelin (1976) 
and the gray by Taylor, Cordes (1993), which is modified 
for the better fit of the kinematic distances by Downes 
et al. (1980). Taylor, Cordes (1993) introduced a "bump" 
structure in Sagittarius arm, but our parallax measure- 
ment does not support it. As Sato et al. (2010a) proposed, 
the "bump" may be caused by the errors of kinematic dis- 
tances. 

It might be possible that MSXDC G034.43±00.24 has 
a large peculiar motion in our Galaxy. When we as- 
sume the flat rotation model, the Galactic rotation Vq = 
220 km s _1 and the distance to the Galactic center Rq = 
8.5 kpc, the radial velocity of this source is calculated 
to be wlsr = 19 km s" 1 . The observed radial veloc- 
ity of the C 18 line is about 58 km s (Rathborne et al. 
2005), which is less than the calculated velocity by about 
40 km s _1 . This might be the line-of-sight component of 
the peculiar motion. Furthermore, comparing with other 



star-forming regions at almost the same galactic latitude 
such as G35.20-0.74 and G35.20-l.74, whose distances 
are measured from parallaxes by Zhang et al. (2009) of 
2 - 19 -o.lo k P c and 3 -27toi2 kpc, respectively, MSXDC 
G034. 43+00. 24 is closer to the Sun. More annual paral- 
lax measurements of Sagittarius-Carina-arm sources with 
VERA reveal the structure of Sagittarius-Carina arm. 

6. Summary 

We carried out the phase-referencing VLBI observa- 
tions of an infrared dark cloud MSXDC G034.43+00.24 
with VERA, and derived its annual parallax of 0.643 ± 
0.049 mas. This corresponds to the distance of 
1.561q \i kpc, which is less than the half of the previ- 
ous kinematic distance of 3.7 kpc. We showed that we are 
able to measure the parallaxes of low-declination source at 
^J20oo ~ +1° with the fitting of right- ascension data only. 
The stable components of the earth's atmosphere cause 
the large scatter in the declination direction of our result. 

Since the distance becomes shorter, the masses of mil- 
limeter cores decrease by one order of magnitude. The 
mass of the most massive millimeter core changes from 
the previous estimations of ~ 1000A/ Q to hundreds M Q . 
The spectral types derived from the luminosity are BI- 
BS, which are still early B types. The position on Galactic 
plane is the closer part of Sagittarius-Carina arm. The ra- 
dial velocity calculated from the flat rotation model has 



Annual Parallax of MSXDC G034.43+00.24 



11 




Fig. 10. Positions on the Galactic plane of some sources VLBI parallaxes are measured, overlaid the the model of the Galaxy by 
Georgelin, Georgelin (1976) (black) and Taylor, Cordes (1993) (gray). "GC" is the position of Galactic center. The open square 
indicates the position of MSXDC G034. 43+00. 24 based on our annual parallax measurements. The position of G14.33— 0.64 (open 
diamond) is calculated from the parallax measured by Sato ct al. (2010a) from the H2O maser observations with VERA. The position 
of G35.20-0.74 (open triangle), G35.20-l.74 (cross) and W51 IRS2 (asterisk) are also based on the parallaxes by Zhang et al. (2009) 
and Xu et al. (2009) with the VLBA 12-GHz methanol maser observations. The position of W51 Main/South (open circle) is based 
on the parallax measurements by Sato et al. (2010b) with H2O maser observations with the VLBA. Uncertainty of distances to W51 
IRS2 is shown by the dotted line. The distance uncertainty for the other sources is too small to be hidden by the symbols. [This 
figure is made from the figure in Sato et al. (2010a).] 



a large offset from observed radial velocities by around 
40 km s , which might be the line-of-sight component of 
the peculiar motion. 

Appendix 1. Details of Data Reduction 

Here, we describe the details of our data reduction. 
Figure 11 shows the procedure of our data reduction and 
the setting parameters in the AIPS software. 

In steps (1) and (2), we read raw data from correlator 
(FITS files) to AIPS. The AIPS task INDXR (Hereafter, the 
CAPITALS OF TYPED FONTS mean the task or verb (pro- 
grams) names of AIPS.) was needed for the preparation 
in the AIPS program, cparm = , , . 02 , denotes 
that the time interval of data calibration is 0.02 minutes. 
[Hereafter, small letters of typed fonts denote the 
adverbs (setting parameters) of AIPS in the format of 



'(adverb names) = (setting value)'.] 

Steps from (3) to (8) were the amplitude calibration 
for all sources. ACC0R (3) makes the normalizations of 
visibilities using auto-correlation data. The intervals and 
integrated durations of this solution were set to 2 minutes 
with solint = 2. After applying this result with CLCAL 
(4), we carried out the amplitude calibration with the sys- 
tem temperatures and antenna efficiencies. We used sys- 
tem temperatures supplied from Mitaka correlator side. 
Some antenna efficiencies for different separation angles 
are supplied from Mitaka correlator, but we used antenna 
efficiencies at the zero separation angle because there is 
little change with separation angles 1 . We deleted all gain- 
curve (GC) tables (inext = GC and invers = -1) with 
EXTDEST (5) and read new one from the text file of the 

1 VERA Status Report, available from 

http : //veraserver . mtk . nao . ac . jp/restricted/index-e . html 
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Target source's data 
(na rrow-band d ata) 

(i) r 



(81 



ID) 



110) 



(11) 



| FITLD | 






| INDXR | 






I ACCOR | 






| CLCAL | 






|EXTDEST| 






I ANTAB | 






| APCAL | 






| CLCAL | 






| BPASS | 






I TBIN | 






| SNEDT | 






| CL( 


:al I 



Reference source's data 
(wide-band data) 



outclass = 'FITLD' 
[+ FITLD CL1] 

inclass = 'FITLD'; cparm = 0, 0, 0.02, 

[+ FITLD SN1] 

inclass = 'FITLD'; solmt = 2 

[+ FITLD CL2] 

inclass = 'FITLD' 

sources — '; interpol — '; calsour =' ' 
snver = 1; gainver = 1; gainuse = 2 

[- FITLD all GCs] 

inclass = 'FITLD'; inext = 'GC; invers =—1 
[+ FITLD GC1] 

inclass = 'FITLD'; tyver = 0; gcver = 1 
[+ FITLD SN2] 

inclass - 'FITLD'; tyver = 1; gcver = 1 
snver = 2; aparm = 0; solint = 

[+ FITLD CL3] 

inclass = 'FITLD' 

sources — '; interpol — '; calsour =' ' 
snver = 2; gainver = 2; gainuse = 3 

[+ FITLD BP1] 

inclass = 'FITLD'; calsour = 'J1800+38' " 
docalib = —1; solint = —1; bpassprm = 1, 0; 
outvers = 1; smooth = 1, 0; 
timerang = (8 min. at the highest elevation) 

[+ FITLD SN3] 
outclass = 'FITLD' 
reading the recalculation file 

[+ FITLD SN4] 

inclass - 'FITLD'; inext - 'SN'; invers - 3 
editting large delay values 

[+ FITLD CL4] 

inclass = 'FITLD' 

souces = "; interpol = 'AMBG'; calsour = " 
snver = 4; gainver = 3; gainuse = 4 



T 



Target source's data 
(narrow-band data) 



(continued from lower left) 



Reference source's data 
(wide-band data) 

| FITTP | (21) 

inclass = 'REF' 



self-calibration and 
imaging with DIFMAP 



[+ REFMAP whole file] | FITLD~| (22) 
outclass = 'REFMAP' 
reading self-calibrated images 



[+ AVER SN4] | FRING~| (23) 

inclass = 'AVER'; calsour ='J1855+02' " 
docalib = 2; gainuse = 2; doband = 
in2class = 'REFMAP; invers = K ncomp = 20, 
flux = 0; nmaps = 1; refant = 1 
solint = 4; solsub = 10; aparm = 0, 0, 0, 0, 1, 2, 3, 
dparm = 1, 10, 30, 0; snver = 4 

[+ AVER SN5] | SNEDT~| (24) 
inclass = 'AVER'; inext = 'SN", invers = 4 
editting large delay data 



(25) | TBIN | [+ FITLD SN5]<- 



frequency conversion 



outclass = 'FITLD' 



mclass = 'AVER' 
inext — SN'; invers = 2 



- | TBOUT 



(26) | TBIN | [+ FITLD SN6]*- 



frequency conversion 



outclass = 'FITLD' 



mclass = 'AVER' 
inext = 'SN'; invers = 5 



(2X1 







CLCAL 


CLCAL 


CLCAL 







outclass = 'FITLD' 

calibration with 'horn-on-dish' method 
[+ FITLD CL5, CL6, CL7] 

inclass = 'FITLD'; sources = "', interpol = 'AMBG'; calsour = " 
snver = 5; gainver = 4; gainuse = 5 
snver = 6; gainver = 5; gainuse = 6 
snver = 7; gainver = 6; gainuse = 7 



mclass = 'FITLD'; mext = 'AN' 

keyword = 'ARRNAM; keystring = 'VLB A' 



[+ AVSPC whole file] | AVSPC~| (13) 
inclass = 'FITLD'; docalib = % gainuse = 4 
doband = 2; bpver = 1; smooth = 

doacor - l; outclass - 'AVSPC 
avoption = 'SUBS'; channel = 4 



(30) | SETJY | [rewriting FITLD SU] 
inclass = 'FITLD' 

sources ='G34.4MMl' 'G34.4MM2' 'G34.4MM3' 'G34.4MM4' ' 
sysvel = 57.3; restfreq = 2.223E10, 5080000 
veltyp = 'LSR'; veldef = 'RADIO'; aparm = 256, 



[+ AVER whole file] | UVAVG~| (14) 
inclass = 'AVSPC'; docalib = 0; doband = —1; doacor = 1 
outclass = 'AVER'; xinc = 1; yinc = 4; zinc = 1 



(31) |~CVEL | [+ CVEL whole file] 

inclass = 'FITLD'; outclass - 'CVEL' 

sources = 'G34.4MM1' 'G34.4MM2' 'G34.4MM3' 'G34.4MM4' ' 
doband = 2; bpver = 1; gainuse = 



[+ AVER SN1] | FRING~| (15) 

inclass = 'AVER'; calsour = 'J1800+38' "', docalib = -1 
doband = —1; refant = 1; solint = 2; solsub = 
aparm = 0, 0, 0, 0, 1, 2, 4, 0; dparm = 0; snver = 1 



(32) | INDXR | 



mclass - 'CVEL' 



[+ AVER SN2] | SNEDT] (16) 

inclass = 'AVER', mext = 'SN'. invers = 1 
editing large delay data 



[+ AVER CL2] | CLCAL~| (17) 

inclass = 'AVER' 
sources - "• interpol = 'AMBG'; calsour = " 
snver = 2; gainver = 1; gainuse = 2 



(33) | IMAGR | [+ image file] 

inclass = 'CVEL'; sources - 'G34.4MM1' " 
docalib = 2; gainuse = 7; doband = —1; 
bchan, echan = (three channels centered at peak) 
cellsize = le— 4, le— 4; imsize = 512, 512 
rashift, decshift = (values in Table 2) 
niter = 10000 



(34) 1 JMFfT~| 



[+ AVER SN3] | FRING~| (18) 

inclass - 'AVER'; calsour = 'J1855+02' " 
docalib = 2; gainuse = 2 
doband = —1; refant = 1; solint = 4; solsub = 
aparm = 0, 0, 0, 0, 1, 2, 3, 
dparm = 1, 10, 30, 0; snver = 3 



blc, trc = (range of maser radiation, one frequency channel) 
ngauss = 1; domax = 1; dopos = 1; dowidth = 1 
niter = 40 



[+ AVER CL3] | CLCAL~| (19) 

inclass = 'AVER' 
sources = "; Interpol = 'AMBG'; calsour = " 
snver = 3; gainver = 2; gainuse = 3 



[+ REF whole file] | SPLIT~| (20) 

inclass = 'AVER' 
sources = 'J1855+02' ". docalib = 2', gainuse = 3 
doband = — I', outclass = 'REF'; aparm = 2, 



(continuing to upper right) 



Fig. 11. Details of the data reduction for an observing epoch. Names in the boxes denote AIPS tasks and verbs. Some important 
adverbs, which are the setting parameters for AIPS tasks and verbs, are shown below the task names in the format of "(adverb 
name) = (adverb values)". Square brackets ([ . . .]s) show the outputs of the tasks or verbs. Plus signs mean creating new AIPS 
tables or files. Minus signs mean deleting AIPS tables. AIPS files are denoted with the CLASS names only because NAME parts 
are used to distinguish observing epochs. 
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zero separation angle with ANTAB (6). The amplitude cal- 
ibration itself was done with APCAL (7). This result was 
also applied with CLCAL (8). 

Bandpass calibration was done with BPASS (9). The 
calibration was applied to all sources by using the data 
of the calibrator source, QSO J1800+3848 (calsour 
= 'J 1800+38') during 8 minutes at the highest ele- 
vation (solint = -1 and timerang = (8 min. at the 
highest elevation)). We used raw data of the cali- 
brator source without any calibrations (docalib = -1). 
bpassprm =1,0 means that this calibration was done 
only for the amplitudes, not for the phases. Ther is lit- 
tle phase variation within and between the IF(s) because 
VERA uses the digital filter bank for frequency fileterings. 
smooth =1,0 denotes that we apply Hanning smooth- 
ing function, whose width is four channels, before band- 
pass calibration. 

Modification of tracking model in Mitaka FX correla- 
tor was done by reading a text file with TBIN (10). This 
text file contains the difference of tracking model between 
Mitaka FX correlator and CALC3/MSOLV, as explained 
in §3. We edited large delay data with SNEDT (11). This 
large delays were introduced because of the bug of the 
program making text files for TBIN. This calibration was 
applied with CLCAL (12). 

Time and frequency integrations were conducted at 
steps (13) and (14). First frequency integration was car- 
ried out with AVSPC. Four frequency channels were inte- 
grated to one frequency channel (avoption = ' SUBS ' and 
channel = 4). This task created new AIPS files, whose 
visibilities are calibrated with the previous calibrations 
(docalib = 2, gainuse = 4, doband = 2 and bpver = 
1). doacor = 1 means that the auto-correlation data are 
also included in the new AIPS files. Time integration was 
conducted with UVAVG. Data of 4 ses were integrated to 
one time step (xinc = 1, yinc = 4 and zinc = 1). 

Calibration of clock parameters was done with the 
global fringe search [FRING, (15)] of the calibrator source 
(calsour = 'J1800+38' ''). The integration time and 
the time intervals of solutions were 2 minutes (solint 
= 2, solsub = 0). The detection threshold was 4 in 
the signal-to-noise ratio (aparm = 0, 0, 0, 0, 1, 2, 
4, 0). After editing large delay data with SNEDT (16), the 
results were applied with CLCAL (17), interpolating 2tt am- 
biguity of phases using rate values (interpol = 'AMBG'). 
Then the global fringe search [FRING, (18)] was conducted 
for the reference source (calsour = 'J1855+02' ''). We 
set the delay and window of 10 ns and 30 mHz, respec- 
tively (dparm = 1, 10, 30, 0). This result was also ap- 
plied with CLCAL (19). After integrating all frequency 
channels of each IF with SPLIT (20) (aparm = 2,0) and 
outputting to FITS file with FITTP (21), the reference 
source was imaged with self-calibration using DIFMAP. 

For phase referencing, we reran the global fringe search 
for the reference source. After loading the image of the 
reference source with FITLD (22), we carried out the global 
fringe search with FRING (23). The image is used for 
the model of the source structure (ncomp = 20, 0, nmaps 
= 1). The integration time and the solution intervals 



were 4 minutes and 0.4 minutes, respectively (solint = 
4, solsub = 10). Editing large delay data with SNEDT 
(24) , this result and the result of calibration of clock pa- 
rameters were transferred to the target source data (25, 
26). Here we needed to convert phases because they de- 
pend on frequencies in <fi = 2itvt, where (ft, v and r are 
phase, frequency and delay. After loading the calibration 
of 'horn-on-dish' method (27), these results were applied 
(28). 

Calibration of Dopplcr effect in VLBI does not include 
the earth's spin. AIPS conducts this calibration only when 
the ARRNAM (array name) keyword in antenna tables 
is VLBA. So we changed it from VERA to VLBA with 
task TABED (29). The parameters of Doppler calibration 
were set with SETJY by rewriting source (SU) tables (30). 
The frequency channel of wlsr = 57.3 km s _ was moved 
to the channel number 256 (sysvel = 57.3, veltyp = 
'LSR', aparm = 256, 0). The rest frequency was set to 
22.23508 GHz (restfreq = 2.223E10, 5080000). The 
calibration itself was conducted with CVEL (31) for the 
target sources (sources = 'G34.4MM1', 'G34.4MM2', 
'G34.4MM3', 'G34.4MM4' "). INDXR (32) needed for 
creating AIPS internal files deleted by CVEL. 

Imaging of target source was carried out with IMAGR 
(33). The imaged source was MM1 (sources = 
'G34.4MM1' "). The size of a pixel was (1 x 10~ 4 mas) x 
(1 x 10~ 4 mas) (cellsize le-4, le-4). The resultant 
images had 512 x 512 pixels (imsize = 512, 512). The 
maximum number of clean components were 10000 (niter 
= 10000). Measurement of positions and fluxes was done 
by fitting elliptical Gaussians with JMFIT (34). The maxi- 
mum step of iteration in the fitting was 40 (niter = 40). 
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